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ABSTRACT 

In this paper we reconsider the problem of magnetic field diffusion in neutron star 
cores. We model the star as consisting of a mixture of neutrons, protons and electrons, 
and allow for particle reactions and binary collisions between species. Our analysis is 
in much the same spirit as that of Goldreich & Reisenegger (1992), and we content 
ourselves with rough estimates of magnetic diffusion timescales, rather than solving 
accurately for some particular field geometry. However, our work improves upon pre- 
vious treatments in one crucial respect: we allow for superfluidity in the neutron star 
matter. We find that the consequent mutual friction force, coupling the neutrons and 
charged particles, together with the suppression of particles collisions and reactions, 
drastically affect the ambipolar magnetic field diffusion timescale. In particular, the 
addition of superfluidity means that it is unlikely that there is ambipolar diffusion 
in magnetar cores on the timescale of the lifetimes of these objects, contradicting an 
assumption often made in the modelling of the flaring activity commonly observed in 
magnetars. Our work suggests that if a decaying magnetic field is indeed the cause 
of magnetar activity, the field evolution is likely to take place outside of the core, 
and might represent Hall/Ohmic diffusion in the stellar crust, or else that a mecha- 
nism other than standard ambipolar diffusion is active, e.g. flux expulsion due to the 
interaction between neutron vortices and magnetic fluxtubes. 
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1 INTRODUCTION 

The structure and evolution of magnetic fields in com- 
pact stars is of considerable interest. In part this is due to 
the large amount of information available: magnetic field 
strengths have been measured, or more often inferred, for 
many neutron stars, allowing astronomers to build up a de- 
tailed picture of the distribution of fields strengths in them. 
Crucially, there seem to be pronounced systematic differ- 
ences in field strength between different classes of neutron 
stars. The young neutron stars and most other 'normal' pul- 
sars have fields of the order 10 12 G, while the older millisec- 
ond pulsars (MSPs) have much lower fields, of order 10 9 G. 
Similarly (relatively) low fields are believed to exist in the 
neutron stars located in accreting low-mass X-ray binary 
systems (LMXBs). In contrast, magnetars, as observed as 
anomalous X-ray pulsars (AXPs) and soft gamma repeaters 
(SGRs), seem to have exceptionally high field strengths, of 
order 10 15 G. 

The way in which magnetic fields slowly evolve in neu- 



tron stars may well be crucial in understanding the relation- 
ship between these classes of objects. There are indications 
that millisecond pulsars acquire their relatively low field 
strength through a period of accelerated field decay while 
accreting in LMXBs, neatly linking the normal, LMXB and 
MSP populations. 

Even more intriguingly, the structure, strength and, 
most of all, the evolution of the magnetic field may provide a 
unique insight into the state of the stellar interior. Field de- 
cay in the hot crust of an LMXB has been invoked to explain 
the reduction in field strength that accretion is believed to 
bring about, giving inf ormation on the conductivity of t he 
crustal matter (see e.g. lCumming. Arras fc Zweibell 12004 0. 
More exotically, magnetic field evolution has been invoked 
to explain the restless behaviour of the magnetars. As was 
argued convincing ly by iThompson fc Duncan! (| 19951 . Il996l . 
1200 if ), the existence of a 10 15 G strength field explains many 
of the observed features of the AXPs and the SGRs. In par- 
ticular, Thompson and Duncan pointed out that the very 
slow rotation rates of the magnetars imply that the energy 
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required to power their emission cannot be supplied by ro- 
tation, but a 10 15 G magnetic field, evolving on the inferred 
lifetime of the magnetars (of order 10 3 — 10 4 years), would 
provide a suitable energy reservoir. 

In order for this model to be self-consistent, there must 
exist a mechanism capable of causing a 10 15 G magnetic 
field to evolve on a timescale of order 10 4 years. The issue 
of magnetic field evolution had be en considered previously 
by iGoldreich fc Reisenegger] i|l992l ) (hereafter GR92), who 
showed that a variety of processes contribute to the field 
evolution in a neutron star. These processes are (i) Ohmic 
decay, in which the current and therefore the field decays 
through the frequent collisions between the different parti- 
cle species; (ii) Hall drift, in which the magnetic field pro- 
duced by the current acts back on the current itself, produc- 
ing a (dissipation free) rearrangement of the magnetic flux; 
(iii) ambipolar drift, in which the charged particles (elec- 
trons and protons) and the magnetic field move relative to 
the neutral fluid (the neutrons). It was the last of these, 
ambipolar drift, that Thompson and Duncan were able to 
invoke to explain the activity of the magnetars. While not 
important for normal field strength stars, they noted that 
the quadratic scaling of the drift velocity with magnetic 
field strength led to field rearrangement timescales of or- 
der magnetar lifetimes for magnetar-like field strengths. In 
contrast, Ohmic decay is unlikely to be im portant in neutron 
star cores (|Bavm. Pethick fc Pineslll969r i. The Hall drift, al- 
though not a dissipative process itself, might cause a re- 
arrangement of the magnetic field on astrophysically rele- 
vant timescales, assuming a magnetar-strong magnetic field 
(GR92). However, as we discuss in this paper, this effect 
is likely to be suppressed in the presence of proton super- 
conductivity. It should be noted, nevertheless, that the pos- 
sibility remains that Hall and Ohmic processes, working in 
concert, might lead to significant field evo l ution in the stellar 
crust; see e.g. iPons. Miralles fc Geppertl (|2009T ). 

It is precisely this issue of magnetic field evolution in 
a neutron star core that we revisit in this paper. Given its 
apparent significance for the magnetar model, we pay par- 
ticular attention to ambipolar diffusion. Our starting point 
is the classic work of Goldreich and Reisenegger (GR92), 
who modelled neutron stars as being composed of a neutron- 
proton-electron plasma, with field evolution being driven by 
binary scattering processes between the species. The anal- 
ysis of GR92 was essentially at the level of estimating the 
key timescales, rather than calculating explicit fully consis- 
tent solutions. Their analysis h a s rece n tly be en extended by 
iHovos. Reisenegger fc Valdivial <|200Sl , I2010T ). in which the 
magnetic evolution problems were solved explicitly in a sim- 
plified geometry, to monitor the validity of the timescales 
estimated by GR92, and to understand the change in geom- 
etry brought about by the diffusive processes. 

Our aim in this paper is somewhat different. Rather 
than obtaining more explicit solutions to the original model, 
we instead will content ourselves with making simple esti- 
mates in the spirit of GR92, but with the addition of the 
crucial ingredient of neutron superfluidity. There is strong 
evidence, both from theory and observation, that mature 
neutron stars are likely to be sufficiently col d that t he nu - 
clear matter is in a superfluid state, see e.g. ISaulsl l|l989h . 
Our aim is to see what qualitative features result from this 
more realistic level of modelling, and, in particular, to see if 



timescales of order of the magnetar ages naturally emerge. 
That the timescales might change significantly is plausi- 
ble, given that superfluidity is known to suppress particle 
collisions and reactions, but introduce a new mechanism 
for coupling the different fluids, namely mutual friction, in 
which the neutron and charged fluids are effectively coupled. 
This coupling might be brought about by the scattering of 
electrons off the vortices pre sent in the rotating superfluid 
ijAlpar. Langer fc Saulsl lT988y It is the changes in the field 
evolution timescales, and their implications for the magne- 
tar model, that we seek to investigate. 

The structure of this paper is as follows. In section[2]we 
set out the form of the equations to be solved, without spe- 
cialising to any particular model for the fluids themselves. 
In section Owe estimate diffusion timescales assuming that 
the neutron star matter is normal (i.e. not superfluid), es- 
sentially recovering part of the analysis of GR92. In section 
2]we repeat the calculation assuming neutron and proton su- 
perfluidity. Finally, in section [5] we summarise our findings 
and discuss their implications for the magnetar model. 



2 BASIC MODEL 

2.1 Ambipolar diffusion model and 
hydromagnetic equilibrium 

Our hydrodynamical model for the neutron star core is a 
magnetised three-fluid system consisting of neutrons, pro- 
tons and electrons (labelled by the index x = {n, p, e}), al- 
lowing for relative motion with velocities v x . Besides the 
proton-electron electromagnetic coupling, our model ac- 
counts for fluid coupling through binary particle scattering 
or vortex mutual friction (depending on whether some fluids 
are in a superfluid state or not). In addition, our model is 
'transfusive' in the sense that it allows for particle conver- 
sion through reactions. 

Naturally, the model also includes several simplifica- 
tions. In particular, in writing down our superfluid hydro- 
dynamical equations of motion, we will neglect the effect of 
entrainment between neutrons and protons, which is present 
when both species are superfluid. This approximation is jus- 
tified in view of the order-of-magnitude precision of our final 
result^]. 

We also need to clarify to what extent we ac- 
count for the effects of the expected type II proton 
superconductivity in the magnetohydrodynamics of the 

outer neutron star core (for a rec ent discussion see 

Glampeda kis. Andersson fc Samuelsson! l|2010l )). A key ef- 
fect we do include in our analysis is the superconducting 
magnetic force appearing in the equations of motion in the 
place of the familiar Lorentz force, due to the presence of the 
quantised proton fluxtubes. An effect we do not include is 
the participation of the fluxtubes in the mechanics of mutual 
friction coupling between the fluids. This is an important 
omission in our model given that the dynamics of proton 

1 Note however that by invoking vortex-electron scattering as 
a possible mutual friction mechanism we are implicitly assuming 
the existence of entrainment on the mesoscopic scale of individual 
vortices, as it is entrainment that causes the protons to magnetise 
the neutron vortices. 
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fluxtubes and their interaction with other stellar compo- 
nents (e.g. the neutron vortices) may play a central role in 
magne tic field evolution (see e.g. iRuderman. Zhu fc Cher] 
(1998)). The magnetic field properties (equilibrium, secular 
evolution) in a 'full' superconducting model will be discussed 
elsewhere. 

The remaining simplifications adopted in our model 
have to do with the particular nature of ambipolar diffu- 
sion. By definition, this is a slow process (i.e. it acts on 
timescales much longer than any dynamical timescale) with 
the star evolving through a series of 'snapshots', each one 
of them representing a hydromagnetic quasi-equilibrium. In 
turn, and given that the magnetic field and particle reactions 
have a small overall impact on the bulk stellar structure, we 
can assume that at any given time the system represents a 
small deviation from a stationary, chemically equilibrated, 
non-magnetic star in uniform rotation. 

At this point one has to be careful with the precise no- 
tion of hydromagnetic equilibrium. In the case of a single 
fluid star, this notion is unambiguous: it refers to a state 
where the fluid motion is that of rigid-body rotation and 
where there is a balance between the magnetic stresses, the 
fluid pressure and the centrifugal force (if rotation is signif- 
icant). This force balance is the basis of all existing calcula- 
tions on hydr o magn et ic equilibria ( for re ce nt examples, see 
lHaskell et ail <|200SD ; ICiolfi et all (|2009l ); lLander fc Jones] 
120091)). 

However, in more realistic multi-component stellar 
models, transfusion (particle reactions) can occur, and the 
nature of the hydromagnetic equilibrium may change due to 
the ability of particle reactions in driving fluid flow and vice- 
versa (this will become apparent by inspection of the mass 
continuity equations ([SJ below). As will be discussed in de- 
tail, these flows themselves create frictional forces, from the 
drag of one species against another or (in the superfluid case) 
from the mechanism known as mutual friction. If the rate of 
particle scattering is high enough, the frictional forces con- 
nected with these diffusive flows might overcome the forces 
generated by pressure gradients and play a dominant role in 
balancing the magnetic stresses. 

We are therefore led to make a distinction between two 
sorts of slowly evolving equilibria: friction- dominated equi- 
librium and transfusion- dominated equilibrium. Roughly 
speaking, the designation is as follows. At high tempera- 
tures, the rate of particle reactions is high, which tends 
to eliminate chemical potential imbalances and the corre- 
sponding pressure gradients, while the rate of particle scat- 
tering is also high, so the drag force induced by the fluid 
flow is strong. In this high temperature regime, the key 
force balance is between magnetic forces and the frictional 
drag forces; hence the designation friction-dominated equi- 
librium. Conversely, at low temperatures, the rate of par- 
ticle reactions is low, so that chemical potential perturba- 
tions and the corresponding pressure gradients persist and 
are therefore important, while the rate of particle scatter- 
ing is also low, so the drag force associated with the diffu- 
sive flow is weak. In this low temperature regime the key 
force balance is between the magnetic forces and the pres- 
sure gradients connected with the chemical potential imbal- 
ances; hence the designation transfusion-dominated equilib- 
rium. In our discussion below we will encounter both types 
of quasi-equilibria. 



Of course, no previous studies have calculated such 
multi-component equilibria. The closest things that exist ar e 
the multifluid but non- magnetic equilibria of iPrbd (|2004h ■ 
who computed stationary equilibria of rotating two-fluid 
stars. The two fluids were both rigidly rotating, and both 
particle reactions and particle scattering were 'switched off', 
so as to allow solutions that were both stationary and did 
not involve frictional forces, only pressure gradients, rota- 
tion and gravity. We will refer to such simplified equilibria 
(magnetised or not) as 'ideal' equilibria. They are of signif- 
icance in as much as they are the sorts of equilibria that 
have appeared previously in the literature, and can be rel- 
atively easily calculated. Construction of an ideal magnetic 
equilibrium would be a useful first step to calculating the 
more physical equilibria. Note that the ideal equilibrium is 
relatively close to the transfusion-dominated equilibrium de- 
scribed above, in the sense that in both, it is pressure forces 
rather than frictional forces that are key in maintaining equi- 
librium. 



2.2 Basic formalism 

We can now assemble the equations that describe our sys- 
tem. According to the model defined above, when we formu- 
late the relevant hydromagnetic equations we can (i) ignore 
time derivatives (effectively filtering out wave phenomena) 
and (ii) work at leading order with respect to the deviation 
from a non-magnetic background. These same approxima- 
tions were also used by GR92. 

First we have the Euler equations, on a per-volume ba- 
sis, referred to a frame rotating with the star. We write down 
one equation to describe the neutrons, and one to describe 
the combined proton-electron charged fluid. We will take 
as our background solution a rigidly rotating non-magnetic 
star. Using standard cylindrical coordinates {zu, z, ip} with 
the angular velocity taking the form Qz, where the hat in- 
dicates a unit vector, we have: 
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where p x are the mass densities (with n x = p x /m x the corre- 
sponding number densities), /i x are the chemical potentials 
and <E> is the gravitational potential. In the proton-electron 
equation we used charge neutrality n p = n c and neglected 
the electron mass (hereafter we use m p = m n = m). 

The perturbed hydromagnetic system is described by 
the pair of Euler equations (adopting the so-called Cowling 
approximation which entails <5<J? = 0), 
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(4) 



where _F cp i and _F mag are, respectively, the coupling and mag- 
netic forces per unit volume. Their particular form depends 
on the state of matter (normal or superfluid) and will be 
specified later. 

The Euler equations are supplemented by the mass con- 
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tinuity equations 

V ■ (n x u x ) = T x (5) 

which feature the particle reaction ratefl F x . Baryon and 
charge conservation dictate that T n = — T p and T p = T c . 
These properties are indeed satisfied by beta decay in neu- 
tron star cores, typically operating through the modified 
Urea reactions. The background system, as described by 
eqns. JTJ) and P). is characterised by chemical equilibrium, 
i.e. p, n = /i p + fi c . This is no longer true in the magnetically 
perturbed system. The departure from chemical equilibrium 
is quantified by the 'chemical imbalance' parameter: 

P = Sfi p + 8p, c - Sp. n , Sjl x = ^2L (5) 

m 

For a small /3 we can write 

T n = mX/3 (7) 

where the coefficient A is a function of density and temper- 
ature. 

Note that for the ideal equilibrium defined earlier, the 
Coriolis terms in the Euler equations would be zero (all fluids 
rotate rigidly), the coupling force F cp i would be set to zero, 
and particle reactions eliminated by setting A = 0. 

Returning to the Euler equations ([3]) and © we can 
combine them and produce an equivalent (but slightly more 
useful) set of equations. Taking the difference between the 
two equations and manipulating slightly gives 

2fi(z x i5) + V/3 + — !— F cpl = — F mag (8) 

XpPn pp 

where x p = p p /p is the proton fraction (p — p n + p p is the 
total density) and 

W = Vp — V n (9) 

We will refer to ([8]) as the 'difference Euler equation'. 
Taking instead the sum of the two equations gives 

29,{Z Xil) + V<5/2 n + XpV/3 = -F ma g (fO) 

p 

where u — v n + x p w is the neutron-proton average velocity, 
weighted by density. We will refer to ()1O0 as the 'total Euler 
equation'. 

By combining the continuity equations in a similar way 
we get the set, 

V ■ {pit) = (f f) 

V ■ (n p w) + n n V (&) ■v B = -(l + m\p (12) 

\n n ) \ n n J 

Strictly speaking, the two terms on the left hand side of 
equation (|f 2p could have comparable magnitudes, unless we 
invoke a weak stratification (in the sense that x p and there- 
fore also n p /n n vary only slightly over length scales of in- 
terest) which would allow us to neglect the second term. 
However, and despite its common use by many authors (in- 
cluding GR92), this is not a well justified approximation in 
realistic neutron star models. Perhaps a better reason for 
neglecting the second term in (|12|l would be to use (|f f 

2 These rates also appear in the Eu ler equation s, in the form of 
'rocket' terms T x v x (see for example iPrbd (l2004t) ) . However these 
are higher order terms in our perturbative scheme. 



which suggests a meridional neutron velocity component 
<p X v n ~ x p (if x w). Following this reasoning we will ne- 
glect the second term in 1)120 . but it should be pointed out 
that none of our (order-of-magnitude precision) results and 
conclusions would change had we kept this term, as long as 
v n < w. Hence, we write 

V ■ (npw) = -A/3 (13) 

where we have defined 

A = (\ + ^ m \. (14) 

Before proceeding to our discussion of ambipolar dif- 
fussion we can make some non-trivial observations based 
on the 'difference' and 'total' equations ((Sj and (1100 in the 
limit of an ideal hydromagnetic equilibrium (i.e. « x = 0, A = 
0, F cp i — 0). Firstly, from eqn. (|8}, we see that a star en- 
dowed with a magnetic field cannot be in chemical equilib- 
rium, as it is the chemical imbalance ft that provides the 
'pressure' gradient V/9 necessary to balance -F m ag/p P , the 
magnetic force per unit mass of charged plasma. At the same 
time, from eqn. 1)100 we see that only in the special case when 
Xp is a function of /3 is -F m a g /p, the magnetic force per unit 
total mass, equal to the gradient of a scalar. This is rem- 
iniscent of the case of a single-fluid barotropic star, where 
the notions of stratification and chemical equilibrium are 
not relevant, an d Fmag/p is always ba lanced by the gradient 
of a scalar (e.g. lHaskell et al.l l|2008h ). Finally, from either 
Euler equation it follows that an axisymmetric system has 
a vanishing azimuthal magnetic component i^ ag = 0. 

2.3 General strategy 

Our strategy will now be to look for solutions to the above 
set of equations, where we regard the magnetic field B and 
therefore also F mag as given, so that the magnetic force on 
the right hand sides of eqns. ((8| and 1)100 play the role of 
source terms, driving departures from chemical equilibrium 
and ambipolar diffusion of the charged particles with respect 
to the neutrons. We will assume that the magnetic field itself 
is axisymmetric, and that the subsequent evolution preserves 
this axisymmetry. This analysis is very much in the spirit of 
GR92-we only seek to make rough estimates of timescales, 
so will not attempt to construct fully self-consistent solu- 
tions, neither for the background nor for the perturbations. 

The departure from chemical equilibrium is given by 
the chemical imbalance ft. We will use the proton-neutron 
velocity difference w as our diagnostic of ambipolar diffu- 
sion, and refer to w as the ambipolar diffusion velocity. Of 
course, ambipolar diffusion refers to the combined motion of 
the charged species (protons plus electrons) relative to the 
neutrons, so for the use of the term 'ambipolar diffusion' to 
be strictly accurate, it is necessary that the electron-proton 
velocity difference w pc = ?T P — v a be smaller than the proton- 
neutron velocity difference: 

w S> Wp C (15) 

We will check if this inequality is satisfied in our calcula- 
tions below. Note that when it is not satisfied, we would 
still expect the charged particle flow to generate magnetic 
field evolution; it is simply that it is no longer appropriate 
to describe it as classical ambipolar diffusion. 
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In the calculations that follow, we will use the pertur- 
bation equations to make rough estimates of the ambipolar 
diffusion velocity w = \w\, which can then be immediately 
converted into timescales (r = L/w, where L is a charac- 
teristic length scale). Within the order-of- magnitude preci- 
sion of this calculation the timescale r is essentially equal 
to the evolution timescale ~ L/v p derived by the magnetic 
induction equation (see, for example, GR92). We will make 
repeated use of eqn. ([8|, the difference Euler equation. As 
noted above, the magnetic force on the right hand of this 
equation is to be regarded as the driving term, sourcing am- 
bipolar diffusion. The terms on the left hand side all depend 
linearly on w. For the Coriolis term this dependence is ex- 
plicit. For the coupling term this dependence is a function 
of the assumed state of matter (normal verses superfluid), 
as will be described below. The term V/3 represents the dif- 
ference in pressure forces on the neutron and charged fluids, 
and can be related to the ambipolar diffusion velocity via 
eqn. (|f 3[) . In the next two sections we explore the ambipolar 
timescales for normal and superfluid neutrons, accounting 
for the relative importance of the different terms as a func- 
tion of temperature. 



3 AMBIPOLAR DIFFUSION: NORMAL 
MATTER 

Before tackling the superfluid problem, we will look at the 
case of normal matter. We will assume that the particle 
reactions are t hose of modified Urea, with an efficiency factor 
l|Sawveril 19891 ) 



A = 2.5 x 10 27 pif T| erg^cm-V 1 



(16) 



where T 8 = T/(10 8 K) and p 14 = p/(10 14 g/cm 3 ). 

The magnetic force can be identified with the usual 
Lorentz force (per unit volume): 



where V = (p/p n )D pn . 

As we will now show, rotation has a minimal impact on 
ambipolar drift motion. Indeed, we can see that the ratio of 
the Coriolis force to the drag force is ~ 20. /T> <§C 1. It turns 
out that the same ratio controls the magnitude of the az- 
imuthal component of the ambipolar drift. To see this, note 
that from eqns. ()2ip and (1221) our assumption of axisymme- 
try implies an azimuthal component F^ag = T>w v + 2Q,w^ 
and F^ag = 2il(uJ 7 + i p w n ), respectively. Equating these 



w v _ 20, ( < 



2Q 

— < 1 23 



which implies a dominantly poloidal ambipolar flow (this 
conclusion assumes that is not exceedingly larger than 
w m , which is a reasonable expectation). Note that if we 
were to consider a non-rotating star then the same argu- 
ment would have led to vanishing azimuthal components 
F v =w t/ '=Q 

mag u - 

We now wish to use (12ip to make an estimate of the am- 
bipolar diffusion velocity. At this point GR92 carried out a 
decomposition on eqn. (I21[) itself, splitting the equation into 
its solenoidal (zero divergence) and irrotational (zero curl) 
parts. This decomposition has the attractive feature that the 
V/3 term is purely irrotational, so the solenoidal projection 
of the equation represents a balance purely between the drag 
and Lorentz forces. However, there is no simple way of de- 
composing the Lorentz force itself into its solenoidal and ir- 
rotational parts, so we will instead make use of a slightly dif- 
ferent decomposition. We will decompose at the level of the 
ambipolar velocity w, splitting it into what we will loosely 
term its 'solenoidal' and 'non-solenoidaP parts 

w = w s + w lls (24) 

where w s and w nB are defined to satisfy 

V ■ {n p w s ) = V • (n p w ns ) = -A/3 (25) 



Fmag = — (V XB)XB 
47T 



(17) 



The coupling forces arise from binary particle collisions, so 
that the force is a function of the particle collision frequen- 
cies £> xv : 



F cp l = PpVpnW + p e T> m We 



(18) 



where w = v p — v n as defined previously, and w cn = 
Ve — Vn ■ These collision frequ encies have been calculated by 
lYakovlev fc Shalvbkovl {l990), who found that the neutron- 
proton scattering (mediated by the strong nuclear force) was 
much more frequent than the neutron-electron scattering 
(mediated by the weak nuclear force), so we will include 
only the former in our coupling force: 



where 



F cp i « ppXViW 



16 rp2 -1/3 



D p „ « 6.6 x 10 1D Tip 



The difference and total Euler equations then become 

2fi(i xw) + Vw + V/3 = — F mag 
Pp 



2Q,(z x u) + VSp n + x p VP = -F mag 



(19) 
(20) 

(21) 
(22) 



With these formal definitions, the non-solenoidal part of the 
velocity field is that part of w which is sensitive to particle 
transfusion, while the solenoidal part is not associated with 
perturbations in the chemical potentials/pressure. When in- 
serted into (|2ip . and dropping the Coriolis force, this gives 



- V 



iv 

A 



/ -» ns \ 

(n p w ) 



= — F n - 
Pp 



(26) 



To make progress we can exploit that fact that both the 
drag coefficient T> and the reaction parameter A « mX are 
temperature dependent, allowing us to compare the relative 
sizes of the drag and pressure gradient forces acting upon 
the non-solenoidal flow: 



Vw ns 
V/3 



D pn \m 2 L 2 
Pp 



5 x 10" 



2/3 
Pl4 ^01 



(27) 



where we used the normalisations Le, — L/(10 6 cm) 
a; p /0.1. This ratio becomes unity at a temperature 



xoi 



a = 1 



4 x 10 



,8 P 



1/12 1/8 



14 



01 



-1/4 



K 



(28) 



Given the steep temperature dependence of the ratio 
q, it makes sense to consider the low and high temperature 
cases separately. If we begin with the T > T a regime, from 
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eqn. (I27|l we see that we can neglect the j3 'pressure' gradient 
term in (126 p to give 



-Fm 



Pp 



(29) 



We see that at high temperatures, the drag forces dominate 
the pressure gradients and ambipolar diffusion represents 
a balance between the drag force and the Lorentz force. 
We therefore see that for T > T a we are in the friction- 
dominated regime of section [2~T1 or, equivalently, the star is 
in a friction-dominated equilibrium. In this regime the de- 
composition (I24|l and (|25[1 is no longer useful and we can 
estimate a common diffusion timescale 



L 4Tvp p L A D pn 

T r^J r-Nj ■ 

w B 2 



3 x 



W a xoipl{ 3 L 2 6 j£-yr. (30) 



where B 1B = B/(10 15 G). 

For T <T a , eqn. (|27p shows we can instead neglect the 
effect of drag on the non-solenoidal flow, so that eqn. (|26p 
becomes 



Vw s - V 



A 



(lipW 1 



— F m 

Pp 



(31) 



We now need to be careful in extracting a diffusion timescale. 
In fact we can calculate two distinct timescales, by following 
GR92 and identifying that part of the Lorentz force that 
represents a departure of -F mag /p p from a perfect gradient: 



Frr. 



PpV/3 



(32) 



Recall that in the single fluid case, the force -F ma g/p was 
purely irrotational. Unfortunately, as discussed in section 
12.21 there does not yet exist a calculation of the magnetic 
field structure of a multifluid star, so we cannot be sure 
how -F mag compares in magnitude with the full Lorentz 
force. For a generic initial stellar field, it is conceivable that 



\Fm& g \. Another possibility (once the system en- 



ters the T < T a regime) is that the flow w s , and therefore 
Jmag, is quenched as a result of short-timescale, dynamical 
evolution; if viable, this process cannot be described by our 
model. Given our ignorance on this point, all we can say is 



that \T n - 



< \Fm 



and we will parameterise accordingly 



in the equations that follow. 

With this in mind, one timescale is associated with the 
non-solenoidal flow maintained by the non-zero (3 which is 
responsible for balancing the magnetic force. This is de- 
scribed by 



A 



Pp 



(33) 



which leads to an ambipolar diffusion timescale 

± ^fL « 6 x 10 9 Xllp ^ yi 

The second timescale is associated with the solenoidal flow 
sourced by ^ ma g, 



(34) 



Pp 



(35) 



Given that |.Fmag| < |-Fmag| we obtain a lower bound on the 
associated diffusion timescale: 

L p p VL 



J~ 11". 



>3 x 10 s Xol p% 3 Ll^-yx 



14 -"0 o2 
15 



(36) 



The timescales given above are in agreement with the ones 
obtained by GR92, aside from our solenoidal timescale (|36|) 
being a lower bound rather than an estimate (GR92 implic- 



itly assumed \ J- m 



F n , 



This level of agreement is not 



surprising given that, so far, our analysis closely followed the 
one in GR92. 

The following scenario then suggests itself: after the 
birth of a neutron star, once all the fast dynamical perturba- 
tions not considered here have died away, the hot (T > T a ) 
star cools, proceeding through a series of friction-dominated 
equilibria as described by eqn. ()29p . When the star cools 
below T a , the evolution can then be dominated by the 
solenoidal flow of eqn. (|35l) . provided that w s is not sup- 
pressed on a short dynamical timescale. Only on timescales 
longer than that of eqn. (|36|l . after which the solenoidal 
flow is likely to die away, will the evolution proceed slowly 
through a series of transfusion-dominated equilibria as given 
by eqn. (|34|) . The current magnetar activity might then be 
the result either of the combined solenoidal/non-solenoidal 
flow of eqn. (|29[) or the solenoidal flow of eqn. (|35[) . In 
either case, it follows that if the initial stellar magnetic 
field does indeed contain a significant /mag component, the 
currently observed magnetars may well be in a regime of 
friction-dominated equilibrium, in which case all previous 
calculations of th eir equilibrium stru ct ure might be ser i- 
ously in error (e . g. iHaskell et all [^OOSMCiolfi et all l|2009h : 
lLander fc Jones! (|2009l )). 

Finally, we can check if this motion represents ambipo- 
lar diffusion, in the sense of eqn. (|15l) . i.e. we need to verify 
that the proton-electron relative speed w pa is much smaller 
than the ambipolar drift velocity w. The speed w pB can be 
easily estimated from the Ampere law to be 



8 x 10" 



Bis 



(37) 



while, for the high temperature flow of eqn. (|29[) . or for the 
low-temperature solenoidal flow of eqn. (|35p . we have 



10" 



Bi 5 



2/3 



x iL 6 Tgp 14 



(38) 



We therefore see that condition (|15l) is satisfied for stars 
with temperatures of order of a few times 10 8 . However, for 
the non-solenoidal mode of eqn. (|33p . we instead have 



6 x 10 



D2 T rriQ 

-11 -PlS-^S-t 8 

2 4/3 

^-oiPw 



(39) 



which is comparable to the velocity estimate of eqn. (I37P . We 
therefore see that in this case the rapid motion of the protons 
with respect to the neutrons implies that this motion can- 
not be interpreted as ambipolar diffusion in the traditional 
sense. However, as described in section [2~l3l we would never- 
theless expect the magnetic field structure to evolve on the 
timescales estimated. The conceptual difficulty associated 
with the violation of the condition (|15p . although poten- 
tially important for ambipolar diffusion in normal matter, 
becomes irrelevant in the presence of proton superconduc- 
tivity (as we discuss in section |4~2)) . 
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4 AMBIPOLAR DIFFUSION: SUPERFLUID 
MATTER 

In the presence of superfluidity it is clear that the preced- 
ing ambipolar diffusion model needs revision. Once a given 
particle species is in a superfluid state, the rates of physical 
processes where these particles participate (e.g. binary colli- 
sions and chemical reactions) are exponentially suppressed. 
Hence, we need to make the following modifications to the 
normal matter model: 



A -)■ TZ sf A 



(40) 



where TZ pn and TZ s f are superfluid suppression factors. These 
are exponentially small numbers for temperatures well below 
the critical temperatures T cn and T cp for neutron and proton 
superfluidity. In the T <C T cn ,T cp regime, and as long as 
the critical temperatures are not too different, these factors 
can be approxi mated by the following fits to the available 
rigorous results dHaensel. Levenfish fc Yakovlevll200o1 . l200ll ; 



iBaiko. Haensel fe Yakovlevll200ll ). 

TZ pn « 0.01 vly p e -(^+fp>, R sf ^3x 10- 4 y 5 p e~ 2y - (41) 

where y n = 1.188 T cn /T and y p = 1.764 T cp /T (the differ- 
ence in the numerical factors reflects the different type of 
pairing expected in the outer core: 1 So pairing for th e pro- 
tons a nd P2 pairing for the neutrons, see, for example, ISaulsi 
(l989h ). Note that in general TZ np ^ TZ si . 

As a result of the combined superfluid suppression in 
particle reactions and collisions, the notion of a friction- 
dominated equilibrium (section [XT} becomes irrelevant; this 
is demonstrated below in detail. Thus, a superfluid star is 
expected to be close to an ideal hydromagnetic equilibrium 
(i.e. in the transfusion-dominated regime) for any tempera- 
ture below the critical ones. 

Besides the suppression described by (|40[) . superfluid- 
ity also adds a new coupling mechanism between the fluids 
through vortex mutual friction. Mutual friction refers to the 
force mediated by the neutron vortices (which are responsi- 
ble for the star's bulk rotation) as a result of their interaction 
with the charged particles and the magnetic field. The stan- 
dard expression for the resulting coupling f orce (per unit 
volume, acting on the neutrons) is give n by (|Hall fc Vinenl 
1 195fj ; lAndersson. Siderv fc Comeril200r3 ). 



F mf = -2p n Q [B'(z x w) + B{(z 



n] 



where 



B = 



n 



B' = 



TZ 



(42) 



(43) 



1 + TZ 2 ' ~ 1 + TZ 2 

The dimensionless coefficient TZ is associated with the effec- 
tive drag force experienced by a moving vortex. For 'stan- 
dard' mutual friction, i.e. scattering of electrons off the 
vortex's magnetic field, the drag value is estimated to be 
jAlpar. Langer fc Sauld Il988l : lAndersson. Siderv fc Corner! 
l2006h 



TZ^4xW~ 



-> B'^B 2 , B^TZ<.1 



(44) 



In this case the coupling between the neutron fluid and the 
charged particles can be considered weak. The opposite limit 
of strong mutual friction corresponds to 



ft > 1 



B' 



TZ 



(45) 



Strong drag could eme rge as a result of interactions between 
vortices and fluxtubes l|Saulslll989l :l Ruderman. Zhu fc Chenl 
ll99Sl : lLmk|[ 2003). The extreme case of perfect vortex pinning 
corresponds to B' — ► 1 and B — > 0. We will not consider the 
case of perfect pinning in this paper, as it would require a 
model which incorporates the detailed dynamics of neutron 
vortices and proton fluxtubes. These effects are beyond the 
scope of this paper and will be addressed elsewhere. 

The total coupling force is now the sum of the drag and 
mutual friction forces: 



cpl 



(46) 



where the drag force takes on the same functional form as 
in the normal fluid case, aside from the presence of the sup- 
pression factor ft pn : 



-fdrag — p P TZ p i\Ty pn W 



(47) 



The drag force is strongly suppressed at low temperatures, 
so that in a sufficiently cold star the mutual friction force 
will dominate the coupling. However, it is possible that for 
higher temperatures (but still below the critical ones) the 
drag force could dominate the mutual friction, despite the 
suppression of particle collisions. We can quantify this argu- 
ment by means of the ratio 



/ = 



FdT 



F n -_ 



TZ pn D p 11 J ' p 

2Q.B 



5 x 10 iJ ft 



14 ,r, XoiTg 



By..;, \ VlS 



where £> max = max[23, B']. Depending on whether / < 1 
or / > 1 we will refer to the 'cold superfluid' and 'hot su- 
perfluid' regimes, respectively. The transition between these 
two regimes occurs at a temperature Tf at which / ~ 1 
and, due to the exponentially varying ft np , it is very rapid. 
We can estimate Tf by considering the simple case where 
r cn = T cp = T c and adopting the weak mutual friction limit 
((44)) . Using a typical value T c = 5 x 10 9 K in J3TJ we find 
that the interval (say) 0.1 < / < 10 corresponds to 



Tf 
Tf 



2.7- 
2.9 



2.9 x 10 5 K, 
3.2 x 10 8 K, 



P 
P 



10 s 
0.1s 



(49) 
(50) 



illustrating the steep variation of / with temperature, and 
the weak dependence of Tf on spin period. These estimates 
for Tf can be slightly reduced if T c is decreased or T cn < 
r c p (as predicted by some pairing models , see discussion in 
lAndersson. Comer fc Glampedakisl (|2005l ')'). 

The transition from 'hot' to 'cold' superfluid hydrody- 
namics is therefore potentially relevant for relatively hot 
neutron stars. The existing state -of-art cooling calculations 
jPons. Miralles fc Gepper?l2009l 'l suggest that magnetar in- 
terior temperatures are expected to takes values of a few 
times 10 8 K, (at the observed spindown ages 10 3 — 10 4 yr). 
Thus, magnetars could be found in either superfluid regime 
and our modelling should be able to account for both pos- 
sibilities. We will therefore look at each case in turn in the 
following two subsections. 

A final important effect of superfluidity (and in par- 
ticular of type II proton superconductivity) is related to 
the magnetic force -F mag . This force is not given by the 
Lorentz force (|17|) but, instead, is related to the smooth- 
averaged self-tension of the proton fluxtube array. The ex- 
plicit form of the superconducting magnetic fo rce is given in 
iGlampedakis. Andersson fc Samuelssonl (|2010T l . For the pur- 
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poses of this paper it is sufficient to use the rough estimate 

H C \B 



F, r 



4ttL 



(51) 



where H c i « 10 15 G is the critical magne tic field associ- 
ated with the self-energy of a single fluxtube jTillev fc Tillevl 
ll990T ). whose density dependence we will neglect. 



Given that F„ 



^*mag we have 



* mag 
-*- mag 



2D.m 2 L 
^3~ ' 



10" 20 TZ st 



1 + 



LlTl (Is 



1/3 
Pl4 Xoi 



<1 



(58) 



where we assumed that the bracket term is of order unity. 
Finally, by combining ()55[) and (|57|l we obtain, for both weak 
and strong mutual friction, 



4.1 The 'cold superfluid' regime 

We first discuss the regime T < Tf where 'cold' superfluid 
hydrodynamics applies, i.e. / < 1, so that the coupling force 
can be approximated as being entirely due to mutual fric- 
tion. In order to account for vortex mutual friction we obvi- 
ously need to retain a non-zero rotation Q in our equations. 
Also, we will continue to assume an axisymmetric system. 

To begin with, it is convenient to decompose the var- 
ious velocities into poloidal and toroidal components using 
standard cylindrical coordinates. For example, 



w = wp + id tp, 



Wp = W VJ + w z 



(52) 



B 



B w 



W ^> W 



(59) 



In other words, superfluid ambipolar diffusion is dominantly 
toroidal. This is in contrast with the properties of ambipolar 
diffusion in a non-superfluid star, see eqn. (|23|l . 

In order to study ambipolar diffusion within the present 
superfluid model we will mainly rely on the components 
(|55|) and (|56|) of the difference Euler equation, as we did in 
the non-superfluid case. Moreover, we will continue to de- 
compose the ambipolar velocity w into solenoidal and non- 
solenoidal components, as in eqn. (1241) . 

Using the property (|59|) we can approximate (|56|) as 



Inserting these in the Euler eqns. © and (|10p . together with 
the force Q42p . we obtain the difference Euler equation 



20 

— B {w™vj + w v <f>) + V/3 



+2Q 1 



B' 



and the total Euler: 



Pp 



(53) 



2 a 



(Vn + XpW^) tp — (v% + x p w v ) 1 



-He p V/3 = -F mag 
P 



(54) 



The toroidal and poloidal parts of the difference Euler equa- 
tion are: 



29, 



B / B'\ 

— w +11 ] IV 

Xrj \ Xt} 



2S1iz7 



B 



& 



= — (55) 
Pp 



+ V/3 = — F mag (56) 
Pp 



Based on these general equations we can arrive at several 
non-trivial results. Firstly, the z component of (|56|l shows 
that the superfluid system does not allow the establishment 
of chemical equilibrium (/3 = 0) unless F mag = 0. However, 
it is almost certain that a generic magnetic field configura- 
tion does not comply with this requirement. Secondly, we 
can show that the magnetic force is nearly poloidal, i.e. 
-Fmag <C J'mag, which is consistent with the proximity of 
the system (see eqn. (|64[1 below) to an ideal axisymmetric 
hydromagnetic equilibrium (which has -F^ag — exactly). 
Indeed, from (1541) and (15611 we have 



vjV st w v + V/3 = — F ma 
Pp 



(60) 



where the effective 'drag' coefficient V s f takes the following 
form in the limits of weak and strong mutual friction: 



-2a 



90 

^strong ^ ial 



(61) 



To make our strategy maximally transparent, we can elimi- 
nate P from eqn. (|60[) using the continuity equation 



r-7 t -*ns 

V ■ (n p Wp 



(62) 



making the decomposition into 
solenoidal velocity terms explicit: 



solenoidal and non- 



&V st (w* + w™) - V 



1 



V/ -* IIS 

■ (n p w P 



— Fmag (63) 
Pp 



This equation is the superfluid analogue of eqn. (|26p . and 
shows that the structure of the ambipolar diffusion problem 
in the superfluid case very closely mirrors that of the normal 
case: we have a driving term, build out of the Lorentz force, 
on the right hand side, while on the left we have terms linear 
in the ambipolar velocity. 

We can now proceed much as in the normal case, de- 
scribed in section [3] The terms on the left hand side of eqn. 
(|63[) are strongly temperature dependent, so we will exam- 
ine the ratio of the non-solenoidal mutual friction coupling 
force and pressure gradient term: 



a B f = 



V/3 



1(T 18 TU 



10" 



B 



rriO 

is 



(64) 



where we used w^ s ~ — XpW^/B for weak mutual friction 
(the strong mutual friction result differs only by a factor 
x p ). Hence a B f 1 for any relevant temperature and so we 
can neglect the mutual friction term linear in w^ s on the left 
hand side of eqn. (|63[l to give 



F£ mg = 2n P (vZ + x v w ,B ) : 



P -f (57) 



mVsfwl 



V ■ (n p wp") 



= —Fmag (65) 
Pp 
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As in the normal case, care must be taken in extract- 
ing diffusion timescales. For a now with a significant non- 
solenoidal component, the pressure gradient term will dom- 
inate, leading to 



-V ■ 



n-pW-p 



A ma 

Pp 



(66) 



We can now estimate the diffusion timescale for the non- 
solenoidal flow, but now it is important to distinguish be- 
tween the toroidal and poloidal parts: as a result of (|59|) . we 
should expect that the shortest diffusion timescale is associ- 
ated with the toroidal component w^ s . Indeed, specialising 
briefly to the weak drag case, the ratio between the toroidal 
and poloidal ambipolar timescales is 



as non-solenoidal flows generate pressure and therefore com- 
positional changes, which take place on a timescale set by 
the slow particle reactions. The solenoidal timescales are 
much shorter in the superfluid case, as the particle coupling 
forces which inhibit such motion are suppressed. This last 
point means that any solenoidal part of the Lorentz foce 
will be rapidly eliminated, leaving the star in a very slowly 
evolving transfusion-dominated equilibrium. This contrasts 
with the normal fluid case, where the (lower bound) on the 
solenoidal timescale (equation ([36])) was of order 10 5 years, 
and was therefore potentially relevant to the currently ob- 
served magnetar population. The implications of the ob- 
tained timescales for magnetar evolution are discussed fur- 
ther in section [5] 



* in 



^ > 1 
B 



(67) 



Then, combining eqns. (|59|) and (|66[) leads to (again assum- 
ing weak mutual friction) 



Xm 2 BH cl L 



4:TvBp 2 x p 

The resulting ambipolar timescale is 



(68) 



: 6 X 1(T 



10" 
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13 



10 



4/3 

4 I H 1B B 1B Ti yf 



10" 



(69) 

17 is a typical 



where His = flci/(10 16 G), and TZ st 
value for the superfluid suppression factor of eqn. (|41l) for T 
close to Tf. The corresponding strong mutual friction result 
differs only by a factor x p , i.e. r" s — >■ i p r™. From eqn. (|59p . 
the corresponding timescale connected with poloidal flow is 
even longer, so we will not write it down. 

Again in analogy with the normal case, we can consider 
the case of a solenoidal flow; for such a flow, eqn. (|65[) be- 



Pp 



(70) 



where, as before, .Fmag/pp denotes that part of the Lorentz 
force -Fmag/pp that represents departure from a perfect gra- 
dient. Again assuming only that |-F m ag| < -Fmag we can 
make a numerical estimate of the corresponding ambipolar 
timescale: 



> 



BHd 



2 x 10" 



ls\ X01P14LI 



(71) 



The timescale connected with the poloidal ambipolar diffu- 
sion can be estimated as: 



^r s > 2 



10" 



B 



H 15 Bib 



(72) 



where we assumed weak mutual friction. As we pointed out 
earlier, the strong drag result is simply a factor x p shorter. 
In either case, the toroidal timescale is independent of B. 

The ambipolar diffusion timescales JB5J, ([711 . (f72j) are 
among the main results of this paper. Compared to the cor- 
responding results of the normal matter model (section [3]), 
these timescales are significantly different. This is not sur- 
prising given the different nature of fluid coupling (mu- 
tual friction) in the superfluid model. Indeed, these results 
can be understood in simple terms as follows. The non- 
solenoidal timescale is much longer in the superfluid case, 



4.2 The 'hot superfluidity' regime 

We now consider the 'hot' regime T > Tf where matter 
is superfluid but particle collisions still provide the main 
coupling mechanism through the force 

F c pl W PpTZnpDnpW (73) 

Thus, we would expect the analysis of section [3] to be di- 
rectly applicable here, after imposing the suppression (|40|) 
and replacing the Lorentz force with the superconducting 
force 1)51 [I. Recall, however, that in the case of normal mat- 
ter, we were able to safely ignore the effects of rotation, as in 
the difference Euler eqn. (|2ip the Coriolis force was found 
to be much smaller than the drag force. However, in the 
present superfluid case the ratio of the Coriolis to drag force 
is increased by a factor TZ pn '- 



2<> 



2 x 10 



1/3 

-16 Pl4 



7? T 2 
'Vpn-i g 



(74) 



Given that the factor TZpn may be very small, it is therefore 
not immediately clear if rotation will be important in this 
'hot superfluidity' regime. Inspection of ()41[) reveals that for 
most of the T > Tf regime the above ratio remains small. 

Then, we can readily modify the results of section [3] To 
begin with, the ratio a and the temperature T a are modified 
as 



a -> K pn .Ks{ a, T a -> (ftp„ft a f )- 1/8 T 



(75) 



From these expressions we can safely conclude that the 'hot' 
superfluid system will be always in the regime T <^LT a where 
both friction-dominated and transfusion-dominated quasi- 
equilibria could be accessible. 

The normal matter ambipolar timescales (|34[) and (|36|) 
are now replaced by 
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(76) 



(77) 



where lZ np ~ 10 23 is a typical value for the collisional sup- 
pression factor for temperatures close to Tf, as estimated 
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using eqn. (141[) . As it was the case in the 'cold' superfluid 
regime, the evolution of the solenoidal and non-solenoidal 
flows in the 'hot' regime are associated, respectively, with a 
very short and very long timescale. Thus, the system is ex- 
pected to be in a transfusion-dominated quasi-equilibrium. 

Finally we check if the condition w 2> w pc is satis- 
fied. It is straightforward to show that this is the case for 
both regimes of superfluidity. For type II proton supercon- 
ductivity the magnetic field appe aring in the Ampere law 
is th e so-called London field &t, (ICartler. Langlois fc Prix! 



l200d ; iGlampedakis. Andersson fc Samuelsson 2010T ). The 
London field it s elf is given by the standard expression (see 
iTillev fc Tiiievl (|l990h for details), 



e 



10" 



(78) 



The intrinsic weakness of the London field implies a rela- 
tive proton-electron velocity w pe many orders of magnitude 
smaller than the numerical estimate p7[) . Related to this 
last point, the enormous weakening of the electric current 
j — en p w pc in the presence of type II proton superconduc- 
tivity also implies a negligible magnetic field evolution due 
to the Hall drift. 



5 DISCUSSION: IMPLICATIONS FOR 
MAGNETARS 

We now discuss the main results of this paper, consisting 
of the ambipolar diffusion timescales (|69p . (|7ip . (|72[) for the 
'cold' superfluid regime and (|76[) , (|77[) for the 'hot' superfluid 
regime, in the context of magnetic field evolution in magne- 
tars. Ambipolar diffusion could provide a viable mechanism 
for magnetic field evolution in the observed magnetar pop- 
ulation provided the associated timescale is comparable to 
the estimated ages ~ 10 3 — 10 4 yr for these objects. 

Let us first examine the non-solenoidal timescales of 
eqns. (f69|) and (f76)) . As a consequence of the superfluid sup- 
pression of particle reactions (represented by the exponen- 
tially small factor 1Z s i) all non-solenoidal timescales become 
extremely long, much longer than the corresponding nor- 
mal matter timescales (|30p . (|34p and the estimated mag- 
netar ages. An interesting exception could be the 'cold' 
superfluid timescale (|69p if the mutual friction coupling 
could be reduced at the level of B ~ 10 — 4 72. B f . However, 
such a very small value of B is unlikely to occur if mu- 
tual friction is weak, given that the rather well understood 
electron scattering by neutron vor tices leads to a much 
stronger coupling B ~ 10~ 5 — 10~ 4 llAlpar, Langer fc Sauls! 



1988; Andersson, Siderv fc Corned l2006h . In principle, one 
could have a very small B in the opposite limit of nearly 
perfect vortex pinning. However, in the absence of input 
from microphysics, B would have to be arbitrarily chosen 
to be of the order ~ 10 -19 in order for our model to 
produce ambipolar timescales of the order of 10 4 yr. Re- 
gardless of this issue of 'fine-tuning', and as we pointed 
out earlier, the regime of vortex pinning would require 
a more sophisticated superfluid model with detailed vor- 
tex/fluxtube dynamics (e.g. the superc o nduct ing model of 
IGlampedakis. Andersson fc Samuelsson! (|2010T )'). 

A similar situation emerges when we examine the 
solenoidal timescales (|71|) . Q72fl and (|77[) . They clearly show 



that any solenoidal component of the Lorentz force will be 
eliminated on a very short dynamical-like timescale, much 
shorter than the estimated ages of the magnetar population. 
This means that the superfluid case is very much simpler 
than the normal fluid case, where the (lower bound) on the 
solenoidal timescale (equation ()36p ) was comparable to the 
estimated magnetar ages, so that, in the normal fluid regime, 
it was not clear if such flows were physically relevant. 

Thus, we conclude that ambipolar diffusion is unlikely 
to be the driving mechanism for magnetic field evolution 
in the observed magnetars unless we invoke mutual fric- 
tion very much weaker than current modelling suggests. This 
does not sit well with the standard magnetar model, where 
magnetic field diffusion is needed to account for the observed 
activity of the SGRs and AXPs. There do, however, ex- 
ist mechanisms alternative to core ambipolar diffusion that 
might supply the required field evolution. One possibility is 
a com bination of Hall and Ohmic pr ocesses in the stellar 
crust (jPons. Miralles fc GeppertJ [20091. More speculatively, 
it is possible that vortex-fluxtube interactions, not included 
in our model, might provide the crucial evolution mecha- 
nism; such a process will be considered elsewhere. Yet an- 
other possibility is that the inner core is composed of exotic 
matter such as hyperons or quark condensates to which our 
model does not apply. Thus the poorly constrained physics 
of the inner core may allow for processes that lead to mag- 
netic evolution on the required timescale. 

It is possible to consider situations where only one of 
the proton and neutron fluids form condensates without re- 
sorting to further calculation. Let us first consider the case 
where the neutrons are superfluid but the protons normal. 
This would apply e.g. if the magnetic field were to exceed the 
'second critical field' above which proton su perconductivity 
is su ppressed, i.e. if B > H c2 ~ 10 16 G (jTillev fc Tillev! 
1990). In such a regime the absence of entrained protons 
around the neutron vortices would suppress the mutual fric- 
tion coupling force |Feibelmanlll97ll ; ISauls. Stein fc Serene! 
1982), so that the hot superfluid regime considered in sec- 
tion 14.21 would be relevant. The lack of proton superflu- 
idity would also have the effect of decreasing the impor- 
tance of the superfluid suppression on the particle colli- 
sion and reaction rates, so the suppression factors in this 
case would lie s omewhere between the values given in eqn . 
(f4TT) and unity dHaensel. Levenfish fc Yakovlevll2000l . l200ll : 
iBaiko. Haensel fc Yakovlevll200ll '). From eqns. (f76)) and ((77)) 
we see this would have the effect of pushing both the 
solenoidal and non-solenoidal diffusions timescales in the di- 
rection of physical interest. 

Now consider the case where the neutrons are nor- 
mal but the protons superco nducting. Such a model formed 
the ba sis of the analysis of lArras. Cumming fc Thompson! 
(2004), who studied the link between magnetic field de- 
cay and thermal emission. In terms of our analysis, the 
lack of neutron superfluidity eliminates the mutual fric- 
tion force, so that we can simply apply the 'hot superflu- 
idity' model of section 14.21 where the mutual friction force 
was not absent but was negligibly small. We can therefore 
make direct use of eqns. (1761) and (1771), with suitable val- 
ues of the suppression factors jBaiko. Haensel fc Yakovlevl 
120011 ; iHaensel. Levenfish fc Yakovlevl bOOlT) . We find that 
the qualitative picture is unchanged, with the non-solenoidal 
timescale remaining much too long to be of physical interest, 
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while the (lower bound) on the solenoidal timescale remains 
extremely short. 

Our model has a number of obvious weaknesses, the 
most serious of which is the neglect of interactions between 
the neutron vortices and the magnetic fluxtubes. The in- 
clusion of such effects would provide a logical extension 
of the analysis of this paper. Our analysis could also be 
improved by considering particular realistic field configu- 
rations. This would allow explicit calculations of how the 
magnetic field evolves, giving more accurate timescales, and 
information on how the geometry of the field evolves in 
time. Ind eed, such calculations have bee n carr i ed ou t re- 
cently by iHovos. Reisenegger fe Valdivial (|2008l . |2010| ) for 
non-superfluid stars. 

Related to this last point, our work illustrates the need 
for the construction of accurate axisymmetric equilibria, to 
act as backgrounds about which to monitor the field evo- 
lution. However, as we have endeavored to point out, the 
construction of such solutions in complicated by the si- 
multaneous existence of magnetic forces and particle reac- 
tions: the magnetic stresses imply the existence of a non- 
zero chemical imbalance, which in turn implies the exis- 
tence of particle reactions, relative motion of the species, 
and particle scattering. Interestingly, in the case of normal 
matter, our analysis in section [3] makes clear that, if the 
-^mag/pp component of the Lorentz force (the part which is 
not the gradient of a scalar) is significant, then the friction- 
dominated equilibrium that results has particle scattering as 
an important component of the force balance. Such forces 
have been absent from all magnetic equilibrium c alcula - 
tions to date, bot h perturbative (e.g . Haskell et al.l |2008); 
ICiolfi et al l (l2009h; iLander fe Jones i|2009l )) and numerical 
(e.g. Braithwaite fc Nordhmdl l |200€ )). and so these normal 
fluid calculations might be significantly in error. Whether 
such a component of the Lorentz force is relevant, or whether 
it is removed on some fast dynamical timescale by processes 
not considered here, is certainly an issue affecting the normal 
fluid case. Clearly, a scheme for producing pseudo-stationary 
equilibria needs to be developed, for both normal matter and 
superfluid/superconducting stars, which would then provide 
a background for more detailed investigations of secular field 
evolution. 
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